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We derive a new asymptotic formula for the Wigner 9j-symbol, in the limit of one small and eight 
large angular momenta, using a novel gauge-invariant factorization for the asymptotic solution of a 
set of coupled wave equations. Our factorization eliminates the geometric phases completely, using 
gauge-invariant non-canonical coordinates, parallel transports of spinors, and quantum rotation 
matrices. Our derivation generalizes to higher 3nj-symbols. We display without proof some new 
asymptotic formulas for the 12j-symbol and the 15j'-symbol in the appendices. This work contributes 
a new asymptotic formula of the Wigner 9j-symbol to the quantum theory of angular momentum, 
and serves as an example of a new general method for deriving asymptotic formulas for 3nj-symbols. 
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I. INTRODUCTION 

This paper contains new asymptotic formulas for the 
9j-, 12 j-, and 15j-symbols when some quantum numbers 
are large and others small. In particular, we derive in 
full detail an asymptotic formula for the 9j-symbol in 
the limit of one small and eight large angular momenta, 
as well as several other, similar results for the 9j- 12 j- 
and 15j-symbols. The main theoretical tool we use is a 
generalization of the Born-Oppenheimer approximation, 
in which the small angular momenta are the fast degrees 
of freedom and the large angular momenta are the slow 
degrees of freedom. 

In standard applications in molecular physics, the 
Born-Oppenheimer approximation couples a set of elec- 
tronic modes by means of a potential energy matrix, in 
which the matrix elements are all functions of position 
and so commute with one another. In our applications 
the modes of the fast degrees of freedom are coupled by 
a matrix of noncommuting operators. The required gen- 
eralization of the Born-Oppenheimer approximation in 
this case was presented by Littlejohn and Flynn [lj and 
applied by those authors to a semiclassical treatment of 
spin-orbit coupling [|. The method involves diagonaliz- 
ing the matrix of operators by means of the Moyal star 
product, in a pertubation expansion in powers of h. The 
results are interesting geometrically, in that there ap- 
pears a bundle whose base space is the classical phase 
space of the slow degrees of freedom, with a fiber that 
is the Hilbert space for the fast degrees of freedom. The 
fiber bundle carries a Berry's connection, whose curva- 
ture contributes to the symplectic form on the base space. 
This paper makes extensive use of the techniques devel- 
oped inRefs. HHI, and assumes a familiarity with them. 
In addition we note the work of Emmrich and Weinstein 
Q, which reinterpreted the transformations of Littlejohn 
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and Flynn as a deformation of the Weyl symbol corre- 
spondence and placed the whole procedure in a more 
mathematical setting. 

In this paper, as in the spin-orbit problem [|, the small 
angular momentum is represented by exact linear alge- 
bra, and the states are represented by multicomponent 
wave-functions. We use Schwinger's model 0] to rep- 
resent the large angular momenta. Some of the issues 
we face in this paper are non-Abelian symmetry groups 
in the presence of the Born-Oppenheimer approximation 
and the problem of constructing gauge-invariant expres- 
sions for wave functions and matrix elements. 

Part of the significance of this work is a new asymptotic 
formula for the Wigner 9j-symbol in terms of a modified 
form of the well-known Ponzano-Regge formula Q for 
the 6j-symbol. The other part is the demonstration of a 
new technique that is generally applicable to the studies 
of higher 3nj-symbols. 

We will briefly review some of the previous works on 
the asymptotics of the 3nj-symbol, when some of the 
angular momenta are small and others are large. Basic 
references for the definitions and properties of the 3n?'- 
symbols include Edmonds [6J , Biedenharn and Louck |7[ , 
Brink and Satchler Q, and Varshalovich [|. Let us use 
the notation (s, I) to classify the various asymptotic lim- 
its of the 3nj-symbols, where s and I denote the number 
of small and large angular momenta, respectively. For 
a given 3nj-symbol, not all possibilities I + s = inj are 
allowed, since the sum of two small angular momenta 
must be small. For the 3j'-symbol, in 1957, Brussaard 
and Tolhoek [lOj used the Stirling's approximation in a 
sum formula for the Clebsch-Gordan coefficient to derive 
the (1, 2) case. In 1968, Ponzano and Regge @ used intu- 
itive methods to guess the (0, 3) case for the 3j-symbol, 
as well as the (0,6) case for the 6j-symbol. in 1960, Ed- 
monds @ found the (1,5) case, which trivially leads to 
the (2,4) case [ill- In 1999, Watson [HIH derived the 
(3, 3) case for the 6j-symbol, as well as the (4, 5) case for 
the 9j-symbol. More recently, Anderson et. al. [fl [l5|] 



2 



extended Watson's result of the 9j-symbol to the (6, 3) 
case where all the small angular momenta are in different 
rows. The other (6, 3) case where all three small angular 
momenta are in the same row is still an open problem. 
Later Haggard and Littlejohn [16[ found the formula for 
(0, 9) case. Our results in this paper will fill in most of the 
gap for the 9j-symbol by deriving asymptotic formulas of 
the (1,8) case and the (2,7) case. 

Now we give a brief overview of our derivation. Af- 
ter we apply the techniques from Q to the 9j-symbol, 
each multicomponent wave-function consists of a spinor 
factor and a factor in the form of a scalar WKB solu- 
tion. The action in the scalar WKB solution is the in- 
tegral of pdx on a Lagrangian manifold. We introduce 
a nearby Lagrangian manifold, the "67 manifold," which 
is just the A- or B- manifold from [17| . This Lagrangian 
manifold appears in the semiclassical analysis of the 6j- 
symbol, and its action is analyzed in [13]. We perturb 
the unknown action around the action on the 6j man- 
ifold, through a two steps process. First we introduce 
an intermediate manifold using the coordinate transfor- 
mation that gives rise to the gauge invariant coordinates 
from 0- By writing the initial gauge-dependent action 
as a line integral, and perturbing the path for the line in- 
tegral from the initial Lagrangian manifold to a path on 
the intermediate manifold, we can write the original line 
integral as a sum of a gauge-invariant line integral and 
a geometric phase. Then we perturb the gauge-invariant 
line integral again by following nearby Hamiltonian flows 
on the intermediate manifold and on the 6j manifold, 
respectively. This second perturbation further splits the 
gauge-invariant line integral into a known action on the 
6j manifold, and a perturbation. When we combine the 
geometric phase and the perturbation with the spinor 
field on the Lagrangian manifold, we obtain a field of 
quantum rotations acting on a fixed spinor. As a re- 
sult, we obtain a factorization of the multicomponent 
wave function into a known scalar WKB factor, namely, 



those for the 6j-symbol from [l7| , and a field of quantum 
rotations of a reference spinor. Although the quantum 
rotations are difficult to handle by themselves, they par- 
tially cancel out when we take an inner product between 
two multicomponent wave functions. In the case of the 
two 6j manifolds, there are two components in the in- 
tersection set, and only one rotation is required for each 
Lagrangian manifold. It turns out that the spinor inner 
product generates a factor of a Wigner d-matrix, and an 
additional relative phase between the two terms of the 
WKB approximation. In the end, we find an asymptotic 
formula for the 9j-symbol in terms of a Wigner d-matrix 
and a modified Ponzano-Regge formula for a 6j-symbol. 

Now we give an outline of the paper. In Sec. HJ we dis- 
play the main result of this paper, namely, a new asymp- 
totic formula for the 9j-symbol where one of the angular 
momenta is small. In Sec. Mil we describe the represen- 
tation we use for the 9j'-symbol, and express it as an inner 
product between two multicomponent wave-functions. In 
Sec. IIV1 we derive a gauge invariant factorization of the 
wave-functions. In Sec. [V] we use this factorization to 
derive the asymptotic formula for the 9j-symbol. Some 
plots of our formula are presented in Sec. I VII The last 
section contains comments and conclusions. The appen- 
dices contain some generalized results for the higher 3nj- 
symbols. 

II. AN ASYMPTOTIC FORMULA FOR THE 
9j-SYMBOL 

We quote the main result of this paper, namely, a new 
asymptotic formula for the 9j-symbol, where one angular 
momentum, j'3 = s, is small compared to the others. We 
have chosen j'3 to be the small quantum number in the 
main formula below, but any other choice can be reduced 
to this case by the symmetries of the 9_7'-symbol. The 
formula is 



Ji J2 312 
s ji J34 
jl3 324 35 



(-1) 



jl+32+ji+j5+2s+V 



A /(2ji3 + l)(2j 3 4 + l)(127ry) 



cos ( (3i + \) i>i + J ~ S7T + M0i + v <i>i J d s v ^{6) . (1) 



Here the indices on the d-matrix are given by fi = 313— ji 
and v = j'34— 3\. They are of the same order of magnitude 
as the small parameter s. The sum in the argument of 
the cosine runs over the six large angular momenta i = 
1,2,4, 5, 12, 24. Out of the eight large angular momenta, 
these are the ones that do not involve the index i — 
3. The geometric quantities V, ipi, <p±, 4>i, and 9 are 
functions of the vector configuration of the tetrahedron in 
Fig. [I] the construction of which we will describe shortly. 
In exactly the same way as they appear in the Ponzano 
Regge formula [5|, V is the volume of the tetrahedron, 



and each ipi is the external dihedral angle at the edge 
Ji, where ip € [0,7r]. The angle 4>i £ [0, ir] is an interior 
dihedral angle of a tetrahedron containing the edge J\, 
but the tetrahedron in question is not the one in Fig. 
[TJ but rather a related, second tetrahedron in Fig. [5] 
Similarly, the angle 04 6 [0, 7r] is the internal dihedral 
angle at J4 in a third tetrahedron in Fig. [3] The angle 
9 e [0, 7r] is not a dihedral angle, but it is the angle 
between Ji and J4, which is adjacent in Fig. [5] and [3] 
Explicitly, the angles <f>\, (f>4, and 9 are defined in terms 
of the vectors as follows: 
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FIG. 1. The tetrahedron constructed from the six edge 
lengths J r , r = 1, 2, 4, 5, 12, 24. 



J 2 




FIG. 3. The vectors configuration from Fig. [TJ are rearranged 
from Ji + J 2 + J 4 + Js = to J 4 + Ji + J 2 + Js = 0. The 
intermediate vector J24 no longer appears. Instead, we have 

Jl4 = Jl + J 4 - 



in Eq. (UJ, so we cannot use the same procedure above 
directly to construct their vector configurations. How- 
ever, the vectors J r , r — 1,2,4,5 in the tetrahedra in 
Fig. [2] and Fig. [3] are related to those in Fig. Q] by par- 
allel translations, and J14 = Ji + J4, so we can find the 
vector configurations in Fig. [5] and Fig. [3] once we have 
those in Fig. [1] Thus, the quantities V,ipi,<f)i,<f)4,d are all 
functions of the six edge lengths J r , r = 1,2, 4, 5, 12, 24. 



III. A REPRESENTATION FOR THE 
9j-SYMBOL 

A. The 9j-symbol as a scalar product 




FIG. 2. The vectors configuration from Fig. [TJ are rearranged 
from Ji + J 2 + J 4 + Js = to Ji + J 4 + J 2 + Js = 0. The 
intermediate vector Ji 2 no longer appears. Instead, we have 

Jl4 = Jl + J4- 



To use Eq. ([TJ), we construct the vector configura- 
tions of the tetrahedron in Fig. [TJ given the six edge 
lengths J r , r — 1, 2, 4, 5, 12, 24, according to the proce- 
dure in the Appendix of |18| . There, a Gram matrix 
G, which consists of dot products between three vec- 
tors at a vertex of the tetrahedron in Fig. [TJ say J 1; 
J 12, — J5, is expressed in terms of the six edge lengths 
J r = j r + 1/2, r = 1,2,4,5,12,24. Then a singular de- 
composition of G gives a vector configuration of the three 
vectors, and the remaining three vectors are constructed 
from the different pairs of the first three vectors. For in- 
stance, the resulting vector configuration for the values 
(Ji, J 2 , J 4 , Js, J12, J24) = (101,123,88,64.5,68.5,92.5) is 
illustrated in Fig. [TJ Notice the tetrahedra from Fig. [5] 
and Fig. [3] contain the edge length J14, which is not re- 
lated to any of the 9 quantum numbers of the 9j-symbol 



We express the 9j-symbol as an inner product of two 
multicomponent wave-functions, using a variation of Ed- 
mond's definition, (6.4.2) in [6j, as follows, 



3i 32 312 

S H 334 
jl3 J24 35 



(5) 



(b\a) 



[(2ji 2 + l)(2j 3 4 + l)(2jl3 + l)(2j 24 + 

Here \a) and \b) are normalized simultaneous eigenstates 
of lists of operators with certain eigenvalues. We will ig- 
nore the phase conventions of \a) and \b) for now, since we 
did not use them to derive our formula. In our notation, 
the two states are 



\b) 



I\ I2 S Z4 ^5 J 13 J 2 4 Jtot 

jl 32 S ji j 5 ji 3 j 2 4 

Jl J2 S 2 J4 J5 J 12 J 34 Jtot 

jl 32 s j 4 js 312 334 



(6) 



(7) 



In the above notation, the large ket lists the operators 
on the top row, where a hat is used to distinguish an 
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operator from its classical symbol function, and the cor- 
responding quantum numbers are listed on the bottom 
row. 

Now we specify the Hilbert space where the two states 
| a) and \b) live. The states belong to a total Hilbert 
space of five angular momenta "Hi <g> 7^2 <S> %4 <S> %5 <8> T^s- 
Each large angular momentum J r , r = 1,2,4,5, acts on 
its own copy of the Schwinger space of two harmonic os- 
cillators, namely, H r = L 2 (R 2 ) [4(. The small angular 
momentum S is represented by the usual 2s + 1 dimen- 
sional representation of SU(2), that is to say, it acts on 
n s = C 2s+1 . 

We now define the lists of operators in Eq. ([6]) and Eq. 
([7]). First we look at the operators I r , r = 1,2, 4, 5, J 2 2 , 
and J 24, which act only on the large angular momentum 
spaces, H r = L 2 (R 2 ), r = 1,2,4,5, each of which can 
be viewed as a space of wave-functions ip(x r i, x r2 ) for 
two harmonic oscillators of unit frequency and mass. Let 
a rfl = (xrfj, + iprn)/ %/2, and = (x rfl - ip r ^)j\f2, fi = 
1,2, be the usual annihilation and creation operators. 
The operators I r and J r i are constructed from a and 
as follows, 



identity matrix. Similarly, the operators S, S 2 should be 
viewed as matrices times the identity operator on all H r . 
The following equations define J 2 3 , J§ 4 , and Jtot- 

(J 2 3 ) a/3 = [h(h + 1) + h 2 s(s + l)]S a p + 2S Ji • Sap, (11) 

(Jidcp = [h{h + 1) + h 2 s(s + l)}S aP + 2hJ 4 • S Q(3 , (12) 

(Jtot)a/3 = (Jl + J2 + J4 + ^5)^0/3 + ftS Q /3. (13) 

B. The multicomponent wave-functions 

We follow the approach used in [l], Q to find the WKB 
form of the multicomponent wave- functions ip* (x) = 
(x, a\a) and i) b a {x) = (x,a\b). Let us focus on tp^(x), 
since the treatment of tp is analogous. Temporarily, we 
will drop the a superscript. 

Let Di, i = 1, . . . , 10 denote the the operators listed in 
the definition of \a) in Eq. ([5]). We seek a unitary opera- 
tor U, such that Di for alH = 1, . . . , 10 are diagonalized 
when conjugated by U, that is, 



Ir 2 ^ ] ^rti^rfi , Jri ^ ^ y ^ruiPi) fif^rv ; (8) 

where i — 1,2,3, and cti are the Pauli matrices. The 
quantum numbers^',,, t — 1,2,4,5, specify the eigenval- 
ues of both I r and J 2 to be j r and j r (j r + l)i respectively. 

The operators J 2 2 and 3 2 4 that define the intermediate 
couplings of the large angular momenta are defined as 
partial sums of the J r , 



Y, tj l t ,{D i ) a pUp v = {k i ),^ (14) 

a/3 

where A;, i = 1, . . . , 10 is a list of diagonal matrix opera- 
tors, and where we have employed the index summation 
notation. Let be the simultaneous eigenfunction for 
the /i th diagonal entries Xi of the operators Aj, then we 
obtain a simultaneous eigenfunction ip a ^ of the original 
list of operators Di from 



J12— J1+J2, J24— J2+J4— (9) 

The quantum numbers j'12 and j'24 specify the eigenvalues 
of the operators J\ 2 and J24, that is, J12O12 + 1) and 
j24(j24 + 1), respectively. See [l9[ for more detail on the 
Schwinger model. 

Now we turn our attention to the operator S 2 that acts 
only on the small angular momentum space C 2s+1 . Let S 
be the vector of dimensionless spin operators represented 
by (2s + 1) x (2s + 1) matrices that satisfy the SU(2) 
commutation relations 

[S it Sj] = i ^ e i]k Sk ■ (10) 
fc 

The Casimir operator, S 2 = s(s + 1), is proportional to 
the identity operator, so its eigenvalue equation is triv- 
ially satisfied. 

The remaining operators J 2 3 , J| 4 , and Jtot act on both 
the T-L r and T-L s . As in the spin-orbit problem Q, these 
are represented by nondiagonal matrices of differential 
operators. Even for the operators I r , J 2 2 , J24, above, 
they should be viewed as differential operators times the 



Here the index [i is called the polarization index, in the 
same sense used in 0. There is no summation over the 
polarization index. Since we are interested in ip a only 
to first order in h, all we need are the zeroth order term 
U of the Weyl symbol matrix of U and the first order 
symbol matrix Aj of Aj. The resulting asymptotic form 
of the wavefunction ip(x), ( (4.6) in IJ), is a product of 
a scalar WKB part Be lS and a spinor part r, that is, 

^ (x) = B{x) e l s ^' h 7-W {x, p) . (16) 

Here the action S(x) and the amplitude B{x) are simul- 
taneous solutions to the Hamilton- Jacobi and the trans- 
port equations, respectively, that are associated with the 
Hamiltonians \[ ■ The spinor r M is the /i th column of 
the matrix U, 

rM(x,p) = U afi (x,p), (17) 
and where p = dS(x)/dx. 
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Now let us apply the above strategy to our problem. 
The Weyl symbols of the operators I r and J„ are I r — 1/2 
and J r i, respectively, where 



Ir 2 ^ y z rfj,Z rfl , 



and where z r 



Jri — 



%r[j,\&i) [AuZrv ? (1$) 



are the 



•£r/i T" VPrii and 2^ r ^ — ^Prfi 

symbols of d and d^, respectively. The symbols of the 
remaining operators have the same expressions as Eq. 
©, (HTJ)-(IT3J), but without the hats. 

Among the operators L>i, Jf 3 and the vector of the 
three operators J to t are non-diagonal. First we look at 
the expression for Jf 3 , Eq. (TTTI) . and see that the zeroth 
order term of the symbol matrix Jf 3 is already propor- 
tional to the identity matrix, so the spinor t must be an 
eigenvector for the first order term Ji ■ S. Let t^(3\) 
be the eigenvector of the matrix Ji • S with eigenvalue 
fill, that is, it satisfies 



Now let us show that all the transformed symbol ma- 
trices of the operators in Eq. ©, namely, the Aj, are 
diagonal to first order. Let us write A[D] to denote the 
operator U^DU, and write A[D] for its Weyl symbol. 
First, consider the operators I r , r = 1, . . . , 4, which are 
proportional to the identity matrix. Using the operator 
identity 



(23) 



we find 



[A(/ r )W = & - V2)<W - {U 0a „ Ir} , (24) 

where we have used the fact that the symbol of a commu- 
tator is a Poisson bracket. Since U atl — t« is a function 
only of Ji, and since the Poisson brackets {3±,I r } = 
vanish for all r = 1,2,4,5, the second term in Eq. (|24|) 
vanishes. We have 



(Ji • S) a /3 



(aO 



(19) 



where /i = — s, . . . , +s. In order to preserve the diagonal 
symbol matrices I r through the unitary transformation, 
we must choose the spinor to depend only on the 
direction of Ji. One possible choice of is the "north 
standard gauge", (see Appendix A of [2]), in which the 
spinor 6 a M is rotated along a great circle from the z-axis 
to the direction of Ji. Explicitly, 



[A(I r )W = (Ir ~ 1/2) V . 

Similarly, we find 

[A( jf 4 )] M „ = j| 4 <5 M „ . 



(25) 



(26) 



Now we find the symbol matrices A(Ji3) for the vector 
of operators J 13, where 



(20) 



where {6i,<j>i) are the spherical coordinates that spec- 
ify the direction of Ji. Note that this is not the only 
choice, since Eq. (flT)|) is invariant under a local U(l) 
gauge transformations. In other words, any other spinor 
t' = e l9< - 3 ^ t that is related to r by a U(l) gauge trans- 
formation satisfies Eq. (Ti~9")) . This local gauge freedom is 
parametrized by the vector potential, 



i( T 00)t 



83 1 



(21) 



which transforms as A^' = A^) — Vj 1 (g) under a lo- 
cal gauge transformation. Moreover, the gradient of the 
spinor can be expressed in terms of the vector potential, 
(see (A.22) in 0), as follows, 



83 1 



-AY 



(aO 



Ji x S 



If 



(22) 



Once we obtain the complete set of spinors /i = 

— s, . . . , s, we can construct the zeroth order symbol ma- 
trix U of the unitary transformation from Eq. (|17[) . 



[A(Jia)U = Ul^iiSafiUpv + hU^S aP Up u . (27) 
After converting the above to Weyl symbols, we find 

[A(Ji 3 )] Ml/ = 3 1 S flv - ihU^iUa^, Ji} + HU^SapUpv 

= Ji^ - itvrM*{T<?\ Ji} + nT^*S a pT^ . 

See (4.10) in Q for a similar calculation. Let us denote 
the ith component of the second term above by T* , and 
use Eq. (|22|) and the orthogonality of r, 



(At)* T ty) 

'a a 



(28) 



to get 



T<=-ihrM*{T<r\j li } 



8tP 



8J 



where in the second equality, we have used the reduced 
Lie- Poisson bracket ( (30) in [19[ ) to evaluate the Poisson 
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bracket {r, Ji}, and where in the third equality, we have 
used Eq. (j2"2")l for dr/dJi. Notice the term involving S 
in T* u cancels out the same term in A(Ji3), leaving us 
with a diagonal symbol matrix, 



[A(Jia)] 



Ji 



See (4.12) in 0. Taking the square, we obtain 



[A(J? 3 )W = (h + »h) 2 6 i 



(29) 



(30) 



Thus we see that A(Jf 3 ) is diagonal. 

Finally, let us look at the last three remaining opera- 
tors J t ot in Eq. (fTBf. Since the symbols J. r for r = 2, 4, 5 
Poisson commutes with Jx, we find A(J r ) = J r . Using 
A(J i3 ) from Eq. (|2"9"|). we obtain 



[A(Jtot)] 
Ji (\ 
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UA { f x J 1 + (J 2 + J 4 + J 5 



(31) 



Therefore all A*, i = 1, . . . , 10 are diagonal. 

For each polarization /i, not counting the eigenvalue 
equation for S 2 , we have 9 Hamilton- Jacobi equations 
associated with the Aj in the 16 dimensional phase space 
C 8 . It turns out that not all of them are functionally 
independent. In particular, the Hamilton- Jacobi equa- 
tions A(Ji) = h ~ 1/2H = jih and A(J 2 3 ) = (h+fih) 2 = 
(ii3 + l/2) 2 h 2 are functionally dependent. The two equa- 
tions are consistent only when we choose the polarization 
(i = ji3 — j\. This reduces the number of independent 
Hamilton- Jacobi equations for S(x) from 9 to 8, half of 
the dimension of the phase space C 8 . These 8 equations 
define the Lagrangian manifold associated with the ac- 
tion S{x) seen in Eq. (|T6|) . 

Now let us restore the index a. We express the multi- 
component wave- function if)® (x) in the form of Eq. (j!6p , 



(32) 



Here the action S a (x) is the solution to the eight 
Hamilton- Jacobi equations associated with the /x th en- 
tries Xf of 8 of the symbol matrices A" . 



h = til 
h = (h 
h = 0*4 
h = (is 



1/2)71, 
l/2)h, 

i/2)h, 
i/2)h, 



(33) 



Jf 4 = (i 24 + i/2) 2 h 2 , 



r(o) 



1 + ^ 
h 



H Ai x J : + (J 2 + J 4 + J 5 ) = , 



and r a = with /i = ji 3 — ji. Note that both r a 
and involve the vector potential Ai, and so are gauge 
dependent quantities. 



We carry out an analogous analysis for ?p b (x). The 
result is 



^ b a (x)=B b (x)e^/ h T b a (x,p), 



(34) 



where 5&(x) is the solution to the following 8 Hamilton- 
Jacobi equations, 



(35) 



h 


= til + 


1/2)71, 


h 


= ih + 


i/2)h, 


h 


= (J4 + 


i/2)h, 


h 


= (is + 


i/2)h, 


J 2 


= (ii2 - 


- \/2) 2 h 2 


T (6) 
J tot 


= J 4 1 


vh 

+ T + 
+4 



h A 4 x J 4 + (Ji + J 2 + Je 



and the spinor r b = t^ v \ v = j 34 — j 4 , satisfies 



The vector potential A 4 is defined by 



r)-r b 



(36) 



(37) 



Again, both r b and Sb involve the vector potential A 4 , 
and so are gauge dependent quantities. 



IV. A GAUGE-INVARIANT FORM OF THE 
MULTICOMPONENT WAVE-FUNCTIONS 

In this section, we derive a gauge-invariant factoriza- 
tion for ip a (x). In general, an action S(x) is defined by 
the integral of pdx on a Lagrangian manifold. The ac- 
tion S a (x) is associated with the Lagrangian manifold 
C a , the level set defined by Eq. (f33|). We will express 
S a (x) in terms of an action on a nearby gauge-invariant 
Lagrangian manifold. First, we pick a reference point zq 
on C a and let xq be the x coordinates of z$ , and write 
S a (x) — f^pdx as a line integral along a path 7 that 
connects xo to x. Because C a is Lagrangian, S a (x) is 
independent of the path 7 chosen. Then we define an 
intermediate nearby gauge-invariant manifold C' a using a 
set of gauge invariant coordinates z' , which we will de- 
fine in the next paragraph. We then express the line 
integral S a (x) on C a as a perturbation of a line inte- 
gral S' a (x) = J^p' (x)dx, where p'(x) is the coordinate 
representation of C' a . Because C a is not Lagrangian, 
S' a (x) will depend on the path 7. Finally, we express the 
line integral S' a (x) on C a as a perturbation of a known 
action S^(x), which is supported on another nearby 
gauge- invariant Lagrangian manifold C 6 J , the so called 
A-manifold from |16| . This manifold appears in the anal- 
ysis of the Ponzano-Regge action of the 6j-symbol. We 
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then show that the perturbations are independent of the 
path 7 chosen, by combining the extra phases generated 
from these perturbations with the spinor field, and show 
that the result becomes a field of quantum rotations of a 
fixed reference spinor on the Lagrangian manifold (x). 
We show the advantages of our factorization by deriving 
a simple expression for the inner products between two 
multicomponent wave-functions. 



A. The intermediate gauge-invariant manifold C! 

The gauge invariant coordinates z' were introduced in 
P, to formulate gauge invariant quantization condi- 
tions, and to show the overall wave function tp(x) in Eq. 
(IT51) is invariant under a local U(l) gauge transforma- 
tion. We use the coordinate map from z to z' here to 
define a gauge invariant intermediate manifold C' a . Note 
the notation z' in [2j denotes vectors of real (x,p) coor- 
dinates, here z' refers to complex coordinates. We define 
the near-identity map Z as follows, 



TjX 



Zrp(Zrp) — Z T p iH[T ) 



\dz, 



rfj. 



(38) 



where Z rfl arc the functions, and z' are the values of 



the functions evaluated at z rll . Since r a only depends on 
Jx, the derivatives dr a /dz rfl for r — 2,4,5 vanish, and 
we have 



z rp Z r p(z r p) Z r p , 



r = 2,4,5. 



(39) 



For r = 1, we use the chain rule, Eq. (JTSJ), and Eq. (|2"Tj) 
to find 



To follow these manipulations, one must pay close at- 
tention to the difference between the functions and the 
values of the functions. 

Looking at the last Hamilton- Jacobi equation from Eq. 
(f33|) . we see that the gauge dependence of J tot is captured 
by the coordinate transformation Z, that is, 



r(«) 



1 



~h 



JioZ+(J 2 +J 4 + J 5 ) = 0. (42) 



Therefore, it is tempting define a new gauge-invariant 
manifold by dropping Z from Eq. (|42[) . This is equivalent 
to treating Z as an active map that maps points on C a 
to a new gauge-invariant manifold C' a that satisfies the 
following Hamilton- Jacobi equations, 



(43) 



h 


- (ii + 


i/2)h, 


h 


- (J2 + 


i/2)h, 


h 


= (j4 + 


i/2)h, 


h 


= (is + 


1/2)H, 


J 24 


= (J24 " 


- l/2) 2 ft 2 , 


(a)' 
tot 


= Jl 1 


+ ^ + (J 2 + J 4 + Js) = o 
h 



It is easy to check that I\o Z = 1\ to first order, and the 
remaining four equations are not affected by Z . 

In other words, if (x,p(x)) G C ai then {x' ,p') = 
Z(x,p(x)) £ L' a and vice versa. Here {x,p{x)) is the co- 
ordinate representation of the Lagrangian manifold C a . 
We let (x,p'(x)) denote the coordinate representation of 
CJ 



B. First perturbation of the action 



dr\ dJi 

Zip, 



z 'ip = Zip(z r p) = z lt , - ifiT^ (Jlj^J 

ihAi ■ ^ & pv zivj 



zip - if 



Using the identity 

Oi aj = Sij + ieijk&k , 
we obtain to first order 

1 



(40) 



Ji 7 - o Z = — z' oiz' 



— A — 

z\ + —A\iZ\ Oi 



1 



1. 



ih 

z\ — —A lk a k Z\ 
1 



-hCiikAuJlk + 7:htijkAikJli 



Jij + h(Ai x Jt) 



(41) 



Now let us relate the line integral S a (x) to the line 
integral S' a (x) — J^p' {x)dx. We express the map Z from 
Eq. (|3"8")l in terms of the x,p coordinates, 



x 'iu = xip + ih{T a ^ 



dx\p ' 

dpip ' 



(44) 
(45) 



To find the value of p'(x) in terms of x, p(x), use the rela- 
tions p' =p'{x') and (x',p'(x')) — Z(x,p(x)), as follows, 



p{x)-ih{T a y 



^L- = p > = p '[x') = p' (x + ih(T a )^ . 



Expanding in powers of h, we find 

Q T a Q T a Qp lv 

dxip dp lv dxip 
dr a 



p[p(x) = Pl p(x) ~ ih(r a y 

= Pl p(x) - ih(r a y 



dx 



(47) 
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and p' rfl {x) — p rf i{x) for r = 2,4, 5. Integrating along the 
path 7, we obtain the relation 

S a (x)=S' a (x)+ih / (r a )t— dxi 
J 1 axi 

= S' a (x) + h [ Ai-dJx, (48) 

where we have used the definition of Ai in Eq. (|21|l in 
the second equality. 

C. Second perturbation of the action 

We now express the line integral S'(x) as a perturba- 
tion of a known action S® 3 (x) on a nearby gauge-invariant 
Lagrangian manifold C^J . The Lagrangian manifold C e J 
is defined by the following Hamilton- Jacobi equations: 

h = (jx + i/2)h, 
i 2 = (j2 + i/2)h, 

h = (H + l/2)h, 

h = (j 5 + i/2)h, 
Jh = Cjm + i/2) 2 ft 2 , 

J tot = Ji + J 2 + J 4 + J 5 = . (49) 

Now we find the difference between S' a (x) and S% 3 (x). 
By lifting 7 to two nearby paths 71 and -f[ on C Q J and C' a , 
respectively, we can model the two nearby paths 71 and 
7J by Hamiltonian flows generated by the Hamiltonians 
in Eq. (|4"9"|) and Eq. (|4"3")h respectively. Since the first five 
Hamiltonians (ii, I2, 14, h, J24) are identical in Eq. (j49|) 
and Eq. (|4Ti|) , we will focus on the segments of 71 and 

7j that are generated by fi • J to t and n • J^ t . Here n is 
some fixed unit vector, and 6 is the parameter along the 
flows of the two segments. Let us denote the projection 
of the two segments onto the x space by x(9) and x'(0), 
respectively. We assume S' a (x'(8)) — S' a (x(9)) holds to 
first order. 

The equations of motion for the Hamiltonian n • J to t 
are 

^ = ^H^ = 2^^ (n " T) -' (50) 

so the line integral is 

S 6j (6) = Im / V z rv dz rv = \ (n- Jtot) d0 = O, (51) 
Jo rv Jo 

where we have used Eq. (|50l) for dz and the definitions 
for J r in Eq. (|18p in the second equality, and where we 
have used J to t = on C^ 3 in the last equality. Thus the 
action generated by part of the flow in 71 vanishes. 



The equations of motion for the Hamiltonian h ■ J[ l 
are the following, 

= - ^ z rp (n ■ a) pv (^1 + _ j - _(n • , 

(52) 

so the line integral is 

S' a (6)= I (A.J^')cW-^ [ (n-3 1 )de 
Jo 1 i Jo 

= I (fi-JiJdff, (53) 
•'1 Jo 

where we have used Eq. (I5"2"j) and Eq. (ITBl in the first 
equality, and the last equation of Eq. (l4U|) in the second 
equality. Thus the part of line integral along 7J does not 
vanish. 

From the difference of the line integrals in Eq. (|5"T]) and 
Eq. (|53|) . we obtain 

S' a {x)=S%(x)-Y,T- I n t -J 1 d9. (54) 

where the sum in the last term is over the segments in 
7J of 71 that are generated by Hamiltonians of the type 
fii • Jtot- 

Combining the two perturbations Eq. ([4"8"]) and Eq. 
([54|) . we can express the action S a {x) on the gauge- 
dependent Lagrangian manifold C a as a perturbation of a 
known action S^(x) on the gauge-invariant Lagrangian 
manifold C & J , 

S a {x) = Sl 3 {x) +hf Ai • dJ 1 - J2 J 1 I *i ■ Ji M ■ 

(55) 

We note that the second perturbation in Eq. ([55)1 van- 
ishes unless 71 contains flows generated by the total angu- 
lar momentum. Moreover, since the only Hamiltonian in 
Eq. (|4^|) that geenerate flows that change the values of Ji 
or J5 is the total angular momentum, the first perturba- 
tion in Eq. (1551) also vanishes unless 7 contains flows gen- 
erated by the total angular momentum. In other words, 
only the flows of the Hamiltonians n • J to t generate the 
two perturbations in Eq. (j55[) . 

Since S a (x) and S® 3 (x) are actions on Lagrangian man- 
ifolds, their values are independent of the path 7 we 
choose to perform this perturbation. Thus, the last two 
terms in Eq. (|55|) together should also be independent of 
the path 71 . We show this is indeed the case by rewriting 
the spinor field on the Lagrangian manifold C & J , which is 
independent of 71 , in terms of these two terms and a field 
of rotational matrices on H j , which is also independent 
of 71, and a fixed spinor. 
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D. Transformation of the spinor field 



We now turn our attention to the spinor T a (z) in the 
multicomponent wave- function from Eq. (|16|) . Let zq be 
the reference point and let 71 be the path from zo to z on 
, as defined above. If 71 is generated by the Hamilto- 
nian flows of (7i, I2, 14, h, ^24); which do not change the 
value of Ji, we have dJi = along 71. Since the spinor 
T a (z) only depend on J 1 , we find 



T a (z)=T a ( Z0 ) 



(56) 



If 71 is generated by the Hamiltonian flows of n • J tot by 
an angle 9, then dJi is nonzero, and is given by 



dJi = n x JxdO . 



(57) 



Now we analyze how a spinor section r a {z), such as the 
one from Eq. (f2"U|). changes as we follow 71. To express 
the gauge dependence of r a (z) on the gauge potential 
Ai, we find it convenient to express r a (z) in terms of the 
parallel transport of the spinor along 71, that is, 



r a (z) = cxp j-i J Ai • dJ x j f a {z 



(58) 



Here f a is the solution to the parallel transport equation 
along 7, 



~d9~ 



J? 



X Ji 



T a {B) 



(59) 



with the initial condition f a (0) = r a (z )- We see that 
the phase in Eq. (1551) cancels the first term in Eq. (f55j). 

To generate a phase that will cancel the second per- 
turbation in Eq. (I55p . we rewrite the parallel transported 
spinor f a in terms of quantum rotation matrices. Let us 
define 



U(n,d0) = e - i( - & - s) 



(60) 



which has the same axis and angle as those in 7. We 
decompose the axis n into components along and per- 
pendicular to Ji, 



n • Ji 

•'1 



Ji 



Ji x (n x 3 X ) 



(61) 



Using the fact that f a (Ji) is an eigenvector of Ji • S from 
Eq. p^|) . and the parallel transport equation from Eq. 
(l59"l) . with dJi defined in Eq. ([5Tj). we find 



E/(n,d<9) f Q (z) 

Ji x (n x Ji) 

n J 



exp < — 1 



■Sd9 



exp \ -i— ji Ji • S d0 I f Q (z) 



exp 



J i 

.Ji x dJi 



7 2 



e l i f a {z + dz). 



where 



o»7 



exp 



• ° ' Jl ^ 



(62) 



(63) 



Moving the exponential to the other side, and integrating 
with respect to d9 along 7 from z Q to z, we find an ex- 
pression for the parallel transport of the spinor in terms 
of the rotation matrix, as follows, 



f a (z) = exp ji/j J dej U(n, 9) T a (z ) . (64) 

In general, if 71 consists of segments generated by overall 
rotations, say 7 = JliTiJ then by applying the above 
formula sequentially, we have 



r'C) = exp <{/> / n(7) / Jl dg 7 [ n^'^^^o) 



Jl 



Putting Eq. (|65j) back into Eq. (|58|). we have 



(65) 



r a (z) = exp <^ / Ai • dJi + i 



■ M 



n(7) • Ji 
Ji 



d9~ 



7i 



r a (z ). 



(66) 



We note that the product of rotation matrices are almost 
uniquely defined by the end points zq and z. To see this, 
project this series of rotations onto a series of classical 
rotations Ri in SO (3). Because both Ji and J5 are in- 
variant under the other Hamiltonians (ii, I2, Ii, I5, J24)) 
the triangle 1-5-24 can only change orientation under the 
flows of the total angular momentum. Thus, this series of 
rotations Ri must take the triangle 1-5-24 at z$ to the tri- 
angle 1-5-24 at z. There is a unique rotation R that does 
this. Thus the product of the rotation matrices f7(rii, 9i) 
is equal to a lift U of the SO(3) rotation matrix R. We 
can write 



U a (x)=Y[U(n l: Bi) 



(67) 



7i 
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We conclude that both the spinor field T a (z) and the 
field of rotation matrices U{z) are independent of the 
path 71. As a result, the phase factors in Eq. (|66|) . and 
the perturbations in Eq. (1551) are independent of the path 
7i- 

When we put Eq. (|55|) and Eq. ([Bo) together 
into the multicomponent wave- function ijj a (x) — 
B{x)e lSa ^/ h T a (x) the phase factors exactly cancel out. 
In the end, we obtain a gauge-invariant representation of 
the wave-function, 

iffi(x) =B a {x)e iS " i ^ h U a {x)T a {x ). (68) 

Similarly, the multicomponent wave-function for the 
state |b) has the following form, 

^(x) = B b (x) e < 3 »/ fi U b (x) t\x ) , (69) 

where for convenience, we have chosen the same reference 
point zq. This is possible, because the reference points 
are arbitrary on both manifolds, and we can choose zq to 
lie in an intersection of C^J and . 



E. the inner products of multicomponent 
wave-functions 

We now take the inner product between the two mul- 
ticomponent wave- functions ip a (x) and ip b ( x )> an d P er ~ 
form a stationary phase approximation. From Eq. (|68[) 
and Eq. (|69)l . we find 

(b\a)=J dxB a (x)B b (x)e i ^ j ^- s " S( - x ^/ h (70) 
(T b (xo)^[U b (x)^U a (x)T a (x ). 

Let Mk denote the k th component of the intersection set 
of the two Lagrangian manifolds C 6 J and C b . Choose 
the reference point z n £ A^o- We know the difference 
of the two actions S b (xk) — S a (xk) is constant in each 
component Mk, since the integrals of pdx cancel in the 
intersections. In particular, since S b (xo) = S a (x ) = 0, 
it vanishes in Mq. Similarly, by their unitary property, 
the product of the rotation matrices cancel out in the 
intersections, and is constant in each Mk- In particular, 
it is the identity matrix in Mq. 

Suppose we know the result of the stationary phase ap- 
proximation to the integral in Eq. (|70l) , when the spinors 
and rotation matrices are absent, is given by 

J dxB a (x)B b (x)e^^ s » 3 ^ h (71) 
= e iK J2 n " ex PWSa>fc) - S b 3 (x k ) - fi k n/2]/h} , 

k 



then by treating the spinors as slowly varying and eval- 
uating them at the stationary points, we get 

(b\a) = e iK Y, ^ exp{i[^'(z fe ) - Sf (z fe ) - /i fe 7r/2]/7i} 

k 

(U° k T b (z Q )y (U° a k T a (z )) . (72) 

In the above formula, the sum is over the components of 
the intersection set, and z k is a point in the fcth com- 
ponent, and where U® k is a rotation matrix determined 
from the rotation that takes zq to a point Zk £ Mk in- 
side C & J , and U® k is similarly defined. The formula Eq. 
(1721 is independent of the choice of Zk, because any other 
choice z' k will multiply both £/° J and by the same ad- 
ditional rotation matrix generated from a path from z k 
to z' k , which cancel out in the product {U® k )^U® k - We 
see that our factorization of the multicomponent wave- 
function leads to expression, Eq. (|72j) . for the inner prod- 
uct between two multicomponent wave-functions. 



V. DERIVATION FOR THE 9J-SYMBOL 
FORMULA 

We derive the main result Eq. ((TJ for the 9j-symbol 
by using the general formula Eq. (|T2"|) . The stationary 
phase approximation without the spinors in Eq. (|71l) was 
calculated in an earlier paper on the 6j-symbol [17] . In 
the following, we summarize the results of that calcula- 
tion, and then evaluate the spinor inner products at the 
stationary phase points to derive the inner product of 
the two multicomponent wave- functions, as in Eq. (|72l) . 
From this inner product, we derive the asymptotic for- 
mula for the 9j-symbol. 



A. The relative action at the two stationary sets 

Let us describe the Lagrangian manifolds and the rel- 
evant paths used to calculate their actions. The two 
Lagrangian manifolds C^J and intersect at two sub- 
manifolds Mq and M\. The relative action S^ J (zx) — 
S b \z\) is calculated along a loop 7 = 71 + 72 + 73, as 
illustrated in Fig. 01 First it goes from a point zo in Mq 
to a point Z\ in M\ along a path 71, which is generated 
by the Hamiltonian flows of J| 4 inside C & J . Then it goes 
back to another point in Mo along a path 72, which is 
generated by the Hamiltonian flows of J\ 2 inside . Fi- 
nally it goes back to the starting point zq along a path 
73 through an overall rotation around — J 5 by an angle 
equal to twice the internal dihedral angle </>5 of the tetra- 
hedron in Fig. Q] We can break this loop up into two 
paths 7 a = 71 and 7 b = —73 — 72 that goes from z$ to z\ 
along the two different Lagrangian manifolds. 

By choosing zq as the reference point for both actions 
S%> and Sf j , we have 
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6j 



£ 6 JI 




M 



Mi 



FIG. 4. Illustration for the intersection of the Lagrangian 
manifolds, and for the paths 71, 72, and 73. 



S 6 J(z )-S 6 b j (z ) = 0. 



(73) 



We quote the results of Eq. (62) in [17[ for the relative 
action at the other stationary phase point. It is given by 



S 6 J (zi) ~ S 6 b 3 ( Zl ) = -2^2 Wi + 2^ , (74) 



where ipi are the external dihedral angles of the tetrahe- 
dron in Fig. [TJ and the sum is over i = 1,2, 4, 5, 12, 24, 
and n is some integer to contributes to an overall phase. 
The amplitudes at both stationary point are given by 



^0. 1 — 



1 



2V12^rV 



(75) 



where V is the volume of the tetrahedron in Fig. [TJ Al- 
together, Eq. ([721 becomes 



denote the direction of the vector J4. From Fig. [5j we 
see that \<p\ is the angle between the (Ji, J24) plane and 
the (Ji,J4) plane, so 4> — — 0i , which is defined in Eq. 
@ • Also see Fig. [3J The inclination angle is simply the 
angle 6 between Ji and J4, so it is the same 9 defined in 
Eq. ©. 



Z 



X 




y 



FIG. 5. The configuration corresponding to the reference 
point zq. 



The gauge choices for the spinors at the reference point 
are arbitrary, and they contribute only an overall phase. 
To be concrete, we choose the spinor t° at zq to be the 
standard eigenvector for S z , that is, 



(b\a) 

ire 

+ (U b T b {z )y (£/ a T Q (z ))exp 



(76) 



where e lK is some overall phase, and U a , Ub are rotation 
matrices associated with 7° and 7 b , respectively. 



B. The reference point zq and the spinor products 

We choose the vector configurations associated with 
zq to have a negative signed volume V , where 6V = Ji • 
( J 2 x J 4 ) . This choice of the starting point is in agreement 
with the choice made in [l7|. In addition, we choose a 
particular orientation for the starting point zq , as follows. 
We put Ji along the z-axis, and Ji x J 2 4 along the y-axis, 
so that J24 and J5 lie inside the xz-plane, as illustrated in 
Fig. [S] Let the the inclination and azimuth angles (6, </>) 



<(zo) 



(77) 



For the spinor r b at Zq, we choose it to be an eigenvector 
of J4 • S in the north standard gauge, 

r b a (z ) = e-« a -^ dUO) . (78) 
Taking the spinor inner product, we obtain 



(r b (z Q ))Hr a (z ))=e^-^ d s (6). 



(79) 



To evaluate the other spinor product, we need to first 
find the rotation matrices U a and Ub, which are generated 
from 7 a and 7b, respectively. The path j a — 71 has no 
flows generated by the total angular momentum, so 



The path 7^ = 
around J5, so 



U a = l. (80) 
-73 — 72 contains only an overall rotation 
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U b = U(J 5 ,2<t> 5 ). 



(81) 



The rotation U b effectively moves J4 to its mirror im- 
age J4 across the rrz-plane, which has the direction given 
by (8, (pi). Thus U b T b (zo) is an eigenvector of 3' 4 ■ S, and 
is up to a phase equal to the eigenvector of J4 • S in the 
north standard gauge. We have 



(82) 



where H± in the phase factor is equal to the area of a 
spherical triangle on a unit sphere. See Fig. [6] 



(b\a) 



cos 



(85) 



#1 



where iii is another holonomy for the Ji vector. 

Because the quantities 4>li 04 j H\, H4 depend only 
on the geometry of the tetrahedron, and are independent 
of fi and is, we conclude that the argument in the cosine 
must be linear in \i and v. Equating the two arguments 
of the cosine in Eq. (|84l) and in Eq. (|55|) . we find the 
linear term to be —(fJ.<fii + ^04 ). We find that Eq. 
and Eq. (|55|) become 




(b\a) = 



cos 



(86) 



— + M01 



Using formulas (10.9.2) and (10.9.6) in @ for the spe- 
cial cases of the 9j-symbol for s = and s — 1, and 
through numerical experimentation, we have determined 
the overall phase to be 




FIG. 6. The holonomy H4 is the area bounded by the curve 
from z to ji along a great circle, and the curve from Ub that 
rotates j'4 to j 4 , and the curve rotating j 4 back to z along 
another great circle. 



Therefore, the spinor inner product at the intersection 
Mi is 

(U b T b (z )y (U a T a (z )) = e ivH * e -*C/"0*i d^(0) . (83) 

Putting this back into Eq. (|76[) . we have an asymptotic 
formula for the inner product, 



(b\a) 



cos 



^2 J ^ - ~r ~ ^ + v ( ~r 



(84) 



Using a different choice of the reference point and 
paths, we can derive an alternative expression for the 
inner product, and eliminate the term H4. Let us choose 
a new reference point zq to correspond to an orientation 
where J4 is along the z-axis, and J 4 x J 12 along the y- 
axis. For the paths, we choose 7° = 73+71 and j b = —72. 
Through essentially the same arguments, we find 



(87) 



where the factor (— l) v is necessary to make the expres- 
sion Eq. ([55]) symmetrical between [i and v, because 
d*{0) = (-l^d^ie). We see that the formula Eq. 
(Qj) satisfies the symmetry of the 9j-symbol obtained from 
a reflection of the j's across the diagonal, which corre- 
sponds to swapping the states \a) and \b). 

Again, through numerical experimentation, we found 
the Maslov index to be 



Hi = l- 4s. 



(88) 



Altogether, we have derived a new asymptotic formula 
Eq. ([T]) for the 9j-symbol, which we repeat here: 



Ji J2 J12 
s ji J34 
J13 J24 h 



V(2j 13 + l)(2j 34 + l)(127rU) 



E 



Ji ipi + -7 - sit + ficpx + vfa d* „(0) . (89) 



VI. 



PLOTS 



We illustrate the accuracy of our formula by plotting 
our approximation Eq. (|89[) against the exact 9j-symbols 
in the classical region for the following values of the j's: 
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3l 32 J 12 
S H 334 
Jl3 J24 J5 



51/2 53/2 28 
1/2 47/2 24 
25 27 j 5 



(90) 



The result is shown in Fig. [7] We see that our approx- 
imations are in good agreement with the exact values 
of the 97-symbol in the classically allowed region, away 
from the caustic. However, it gets progressively worse as 
we come closer to the caustic, where V — and the am- 
plitude blows up. We can see this clearly from the error 
plot in part (a) of Fig. [5] 



o 

E 

P. 




FIG. 7. Comparison of the exact 9j-symbol (vertical sticks 
and dots) and the asymptotic formula Eq. (|89|) in the classi- 
cally allowed region, for the values of j's shown in Eq. (|90p . 
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FIG. 9. Comparison of the exact 9j-symbol (vertical sticks 
and dots) and the asymptotic formula Eq. ()89[) in the classi- 
cally allowed region away from the caustics, for the values of 
jr's shown in Eq. (l9lj) . 
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FIG. 10. Comparison of the exact 9j-symbol (vertical sticks 
and dots) and the asymptotic formula Eq. (|89|l near a caustic, 
for the values of j's shown in Eq. ()91[) . 



FIG. 8. Absolute value of the error of the asymptotic formula 
Eq. ((89)) for (a) the case in Eq. (|90)) . and (b) the case in 
Eq. (|9ip . The error is defined as the difference between the 
approximate value and the exact value. 

Since the asymptotic formula Eq. (f59"| should become 
more accurate as the values of the j's get larger, we plot 
the formula against the exact 9j for another example Eq. 
(f9"Tj) . which is roughly Eq. (|90| rescaled by a factor of 4. 
The result in the classically allowed region away from the 
caustic is shown in Fig. ® and the result in regions close 
to the caustic is shown in Fig. [TUJ The errors for this 



case is shown in part (b) of Fig. El We can see that the 
error also scales with the values of the j's. 



31 32 312 ) { 201/2 205/2 89 ) 
s 34 334 } = { 3/2 197/2 99 \ . (91) 
3i3 324 k J I 100 92 j s J 
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VII. COMMENTS AND CONCLUSIONS 



We have computed the asymptotic formula for the 
Wigner 9j-symbol when some angular momenta are large 
while others are small. In so doing, we have formulated a 
gauge-invariant form of the multicomponent WKB wave- 
functions that appears in the definition of the 9j-symbol. 
Our gauge-invariant form retains the factorization of the 
wavefunction into a scalar WKB part and a spinor part. 
Currently, there is no general theory of finding gauge- 
invariant multicomponent WKB wave-functions, so our 
formulation here provides an interesting special case for 
investigating when such a gauge-invariant form is possi- 
ble. 

There are several directions to extend the present work. 
A natural extension is to find the analytic continuation 
of our formula to the classically forbidden region where 
the Lagrangian manifolds do not have real intersections. 
While the analytic continuation of the cosine term in the 
6j-symbol into its classically forbidden region is known, 
it is unclear what we should do with the d-matrix and 
the complex argument 9 in the forbidden region. This 
will be the subject of future research. 

Another extension is to apply our technique to the 12j- 
symbol with two small angular momenta, and to the 15j- 
symbol with three small angular momenta. This is easy, 
because the same Lagrangian manifolds in this paper, 
namely, those that appear in the semiclassical analysis 
of the 6j-symbol in the body of this paper, also appear 
in those cases for the 12j- and 15j-symbols. Since the 
derivations are almost identical, we will simply display 
these results for the 12j-symbol in Appendix [B] and re- 
sults for the 15j-symbol in Appendix [Cl 

It is more difficult to apply our method to the case 
of the 12j-symbol with one small angular momentum, 
or to the 15j-symbol with two small angular momenta. 
The relevant Lagrangian manifolds are those that appear 
in the semiclassical analysis of the 9j-symbol. We will 
present our results for those cases in a future article. 

One limitation of our method is that it is not applicable 
to cases where there are three small angular momenta in 
a single row or column of the 9j-symbol. The asymptotic 
formula for the (6, 3) case when the three small angular 
momenta are placed in the same row is still unknown. In 
such cases, the matrix operators that define the quantum 
states are more complicated. We expect a solution will 
involve a combination of the method developed in this 
paper and the summation formula used in ^2Hl5| . This 
interplay of two very different techniques will be inter- 
esting for a future research project. 



3l S 2 J 12 
S3 U 3U 
313 324 35 



(Al) 



(b\c 



where 



[(2 J12 + l)(2j 3 4 + l)(2ji3 + l)(2j M + l)]i 



ii s 3 js j 13 J 24 Jtot 

jl s 2 S3 j 4 35 313 324 



\b) = 



I\ S2 S3 Z4 I§ 

jl 32 S3 j 4 j 5 ji2 J34 



T2 t 

12 J 34 J tot 





(A2) 



(A3) 



We choose the four spinors just as in the 9j case with one 
small angular momentum. Then the Hamilton- Jacobi 



and Jf 4 , 



equations for the operators J13, J24, J (2-. 
ply pick out the polarizations. The remaining Hamilton- 
Jacobi equations for both states after the perturbations 
are performed become 



Ji = (ji + l/2)S, 
J 4 = Cj4 + l/2)», 
J5 = (j 6 + l/2)ft, 
= Ji + J 4 + J 5 . 

These equations describe a Lagrangian manifold that 
corresponds to a triangle having the three edge lengths 
Jl, J4, J5, illustrated in Fig. [TT] 




FIG. 11. The angle 9 is defined as the exterior angle between 
the edges Ji and J4 in a triangle having the three edge lengths 
Ji, J 4 , ^5- 



Appendix A: The 9j-symbol with two small angular 
momenta 

We can also take the two angular momenta to be small. 
Let us write down the definition of the 9j-symbol when 
S2 = 32 and S3 = js are small. It is 



Because the two Lagrangian manifolds are identical, 
the inner product of the scalar part of the WKB wave- 
functions is unity, and only the spinor products con- 
tribute. After putting in the correct overall phase, we 
get 
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3l S 2 J12 

S3 h hi 

3l3 324 is 

^_l^(j4.+js)+S2+S3+jl2+jl3 

y/(2j 13 + l)(2j 24 + l)(2ji 2 + l)(2j 34 + 1) 



(A4) 



We plot the formula Eq. (|A4I) below against exact val- 
ues of the 9j-symbol in Fig. [121 f° r the following values 
of the j's: 

ji s 2 j 12 ) ( 67 1/2 135/2 
s 3 k 334 } = { 3/2 54 111/2 } . (A6) 
313 324 35 J [ 135/2 107/2 j 5 

Appendix B: The 12j-symbol with two small angular 
momenta 



where 9 is the exterior angle between the vectors J\ and 
J4 in the triangle formed by J\ , J4 , and J5 , as illustrated 
in Fig. [IT] The angle 6 is given by 



1 ( Ji • J4 



Ji J4 



1 I J i + J 4 J ; 
2 Ji J4 



(A5) 




J 5 



FIG. 12. Comparison of exact 9j-symbol (vertical sticks) and 
the asymptotic formula Eq. (|A4[) . for the values of j's shown 
in Eq. $K$ . 



We now treat the case where two angular momenta 
are small. We take j± — s\ and j'5 = S5 to be small. 
Using (A4) in [2(3] for the definition for the 12j-symbol 
of the first kind, the 12j-symbol can again be written as 
a scalar product of two multicomponent wave-functions, 
as follows: 



si 32 312 3125 
33 34 334 3135 

313 324 s 5 j 6 



(b\a) 



{[312) [334} [j'13] [j'24] [7125] [j'135]}' 



where 



\a) 



\b) 



I2 13 14 S§ Iq J12 J34 J125 Jtot 
Si 32 33 34 s 5 J6 312 334 J125 



Si I2 I3 h S| Iq J13 J24 J135 Jtot 
Si J2 33 34 S 5 36 313 324 J 135 



(Bl) 



(B2) 



(B3) 



Using the formalism developed in this paper, we can now 
find a new asymptotic formula for the 12j-symbol up to 
an overall phase. By expressing the 12j-symbol in terms 
of the 9j-symbol in the special cases si — or 55 — 
from (Al) in [2(j, and by using Eq. ([lj, we can determine 
the overall phase. The result is an asymptotic formula 
for the 12j-symbol when two of the angular momenta are 
small: 



Sl 32 312 J125 
33 34 334 J135 

jl3 324 S 5 j 6 



= (-1) 



324 +334 +J125 + 3135 +31 +35 +V\+V5 



V[3l2}[3l25}[jl3}[jl35mTrV) 

(si + s 5 )n + (hi + n 5 )4>2 + {v\ + V5)4>3 



(B4) 



Here m = j 12 - 32, v x = j 13 - j 3 , fi 5 = j 125 - j 12 , 
and 1/5 = j'135 — j'13. The sum in the argument of 
the cosine runs over the six large angular momenta 
i = 2,3,5,6,24,34. The geometric quantities V, ipi, 4> 2 , 
03 , and 9 are related to the tetrahedron in Fig. [TH which 



has the six edge lengths J,; = ji + 1/2, i = 2, 3, 5, 6, 24, 34. 
As before, V is the volume, and each is the external 
dihedral angle at the edge Ji. The angle 4> 2 is the an- 
gle between the plane spanned by (J 2 , J3) and the plane 
spanned by (J2, Je)- The angle 4> 3 is the angle between 
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FIG. 13. The volume V and the external dihedral angles 
ipi are defined on the tetrahedron with the six edge lengths 
J3, J5, J2, Je, J24, J34- 



10" 9 




b 1 1 1 I 1 1 1 1 I 1 1 1 1 L 



65 70 75 

fa 

FIG. 14. Comparison of the exact 12j-symbol (vertical sticks 
and dots) and the asymptotic formula Eq. (|B4[l , in the classi- 
cally allowed region away from the caustics. The values used 
are those in Eq. ()B8|I . 



We illustrate the accuracy of the approximation Eq. 
(|B4[) by plotting it against the exact 12j-symbol in Fig. 
[14] for the following values of the j's: 



si 32 312 3125 ) { 1/2 201/2 100 101 ) 
33 34 334 has > = < 213/2 199/2 117 105 \ . 
313 hi s 5 j e J [ 106 98 1 j 6 J 

(B8) 

We see that the agreement is excellent. 



Appendix C: The 15j-symbol with three small 
angular momenta 

We now treat the case where three quantum numbers 
are small. We take j'3 = S3, j'5 = S5, and j$ — sq to be 
small. By using the spin network definition on page 66 
in [2l[ for the definition for the 15j-symbol of the first 
kind, we write the 15j-symbol as a scalar product of two 
multicomponent wave-functions, 



!3l 32 312 jl25 jl256 I 
S3 34 334 .7*135 jl356 > (CI) 
313 324 S 5 s e j 7 J 

m rj 

{ [il2] [334} [313} [324} [3125} [3135} [?1256] [jl356] } 5 ' 

where 



the plane defined by (J2, J3) and the plane denned by 
(JsjJe). The angle 9 is the angle between J 2 and J3. 
The explicit expression for the angles fa, 03, and 9 are 
given by the following equations: 



(J 2 x J 3 ) • (3 2 x J e 



|J 2 x J 3 | 
(Ja x J 2 ) 



I J2 x J 6 | 
(J3 x J 6 ; 



J.3 X J 2 J3 X Jfi 



J2 • J 3 

J2J3 



(B5) 
(B6) 
(B7) 



h h s§ I4 s 5 



2 S|/ 7 J 



2 t2 t2 
12 J 34 J 125 



256 J tot 



31 32 S3 34, S5 S 6 37 312 334 3125 J1256 







(C2) 



\b) 



h h S 3 14 Si 



S 2 



h 3 



T 2 

J 24 



T 2 T 2 

J 135 J l 



356 J tot 



31 32 S3 31 S 5 S 6 J 7 J 13 224 J135 J1356 







(C3) 

Following the strategy in this paper, we derive an asymp- 
totic formula for the 15j-symbol when three angular mo- 
menta are small, up to an overall phase. The result is 
the following formula: 



31 32 312 J125 .71256 
S3 34 J*34 3135 .7*1356 
jl3 324 S 5 S 6 37 



^ _ ^ \jl +32 +34 +h +2s 3 +U 3 +/J5 +M6 



COS 



\/[7*34] [713] b'l35] [71356] [to] b*1256] (127rV) 

+ K30i - (M5 + M6)^12 - (f5 + f6*)<£l 



(C4) 
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Here fi 3 = j 34 - ji, v z = j 13 - j u fi 5 = ji 25 - ji2, 
»b = Ji35 - 3 13, Me = Ji256 - J125, and v e = ji 356 - J135. 
The angles <f>i and 4>\2 are internal dihedral angles in 
Fig. [TS1 in other words, 4>\ — ir — ipi and <pi2 — tt — ^12- 
Similar to the previous cases, the angle c/4 is the angle 
between the (Ji, J4) plane and the (Ji, J24) plane. The 
angle <f>\ is the angle between the (Ji, J4) plane and the 
(J4, J12) plane. Here we put primes on these angles to 
distinguish them from the internal dihedral angles <j)\ and 
4 . The angle Q\ is the angle between Ji and J4. The 
angle 9 2 is the angle between Ji and 3 12. Explicitly, the 



angles ■ 



and 9 are given by the following equations: 



cos 



cos 



(Jl X J 4 ) ■ (Jj X J 7 ) 
|Jl x J 4 | I Ji X J 7 | 

(J 4 X Ji) ■ (J 4 X J 7 ) 
IJ4 x Ji I |J 4 X J 7 | 



COS 



Jj_jj 

J\Ji 



(C5) 
(C6) 
(C7) 




FIG. 15. The volume V and the external dihedral angles 
ipi are defined on the tetrahedron with the six edge lengths 

Jl, Jl, Ji, Jl, Jl2, JlA- 



We plot the exact values of the 15j-symbol against our 
approximation Eq. (|CJ4|) for the following values of the 
fs: 



Jl ]2 J12 J125 J1256 
S3 ji J34 jl35 jl356 
jl3 hi S 5 S 6 j 7 

203/2 207/2 96 97 98 
3/2 199/2 100 100 101 
101 108 1 1 j 7 

We see that there are generally good agreements. 



(C8) 




FIG. 16. Comparison of the exact 15j-symbol (vertical sticks 
and dots) and the asymptotic formula Eq. (|C4[I . for the values 
of j's shown in Eq. (|C8f) . 



[1] R. G. Littlejohn and W. G. Flynn, Phys. Rev. A 44, 5239 

(1991) . 

[2] R. G. Littlejohn and W. G. Flynn, Phys. Rev. A 45, 7697 

(1992) . 

[3] C. Emmrich and A. Weinstein, Commun. Math. Phys. 
176, 701 (1996). 

[4] J. Schwinger, in Quantum Theory of Angular Momentum, 
edited by L. C. Biedenharn and H. van Dam (Academic 
Press, New York, 1965). 

[5] G. Ponzano and T. Regge, in Spectroscopy and Group 
Theoretical Methods in Physics, edited by F. Bloch 
(North-Holland, Amsterdam, 1968). 

[6] A. R. Edmonds, Angular Momentum in Quantum Me- 
chanics (Princeton University Press, Pinceton, 1960). 

[7] L. C. Biedenharn and J. D. Louck, Angular Momentum 
in Quantum Physics (Addison- Wesley, Reading, Mas- 
sachusetts, 1981). 

[8] D. M. Brink and G. R. Satchler, Angular Momentum 
(Clarendon Press, Oxford, 1968). 



[9] D. A. Varshalovich, A. N. Moskalev, and V. K. Kher- 
sonskii, Quantum Theory of Angular Momentum (World 
Scientific, Singapore, 1981). 

[10] P. Brussaard and H. A. Tolhoek, Physica 23, 955 (1957). 

[11] A. F. Nikiforov, S. K. Suslov, and V. B. Uvarov, Polyno- 
mials of a Discrete Variable ( Springer- Verlag, New York, 
1991). 

[12] J. K. G. Watson, J. Phys. A. 32, 6901 (1999). 

[13] J. K. G. Watson, in The Role of Rydberg States in 
Spectroscopy and Photochemistry: Low and High ryd- 
berg States, Understanding Chemical Reactivity Series, 
Vol. 20, edited by C. Sandorfy (Dordrecht, Kluwer, 1999) 
p. 293. 

[14] R. W. Anderson, V. Aquilanti, and C. da S Ferreira, J. 
Chem. Phys. 129, 161101 (2008). 

[15] R. W. Anderson, V. Aquilanti, and A. Marzuoli, J. Phys. 
Chem. A 113, 15106 (2009). 

[16] H. M. Haggard and R. G. Littlejohn, Classical and Quan- 
tum Gravity 27, 135010 (2010). 



18 



[17] V.Aquilanti, H. M. Haggard, A. Hedeman, N. Jeevanjee, 
R. G. Littlejohn, and L. Yu, e-print arXiv:1009.2811vl 
[math-ph] (2010). 

[18] R. G. Littlejohn and L. Yu, J. Phys. Chem. A 113, 14904 
(2009). 



[19] V. Aquilanti, H. M. Haggard, R. G. Littlejohn, and L. Yu, 
J. Phys. A. 40, 5637 (2007). 

[20] H. A. Jahn and J. Hope, Phys. Rev. 93, 318 (1954). 

[21] A. P. Yutsis, I. B. Levinson, and V. V. Vanagas, Mathe- 
matical Apparatus of the Theory of Angular Momentum 
(Israel Program for Scientific Translations, Jerusalem, 
1962). 



